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ABSTRACT 

Continuing  our  program  to  understand  the  geometry  and  dynamics  of 
floating  four-bar  linkages,  we  explore  the  relative  equilibria  of  an  assembly 
that  admits  symmetric  configurations.  We  show  that  a  symmetric  config¬ 
uration  is  a  relative  equilibrium.  As  we  vary  certain  kinematic  parameters 
which  preserve  the  symmetry,  a  symmetric  relative  equilibrium  is  bifurcated. 

The  type  of  bifurcations  can  be  either  supercritical  or  subcritical  pitchfork. 

The  stability  of  the  relative  equilibria  at  symmetric  configurations  is  in¬ 
vestigated.  Elementary  techniques  of  singularity  theory  are  applied  in  the 
analysis  of  the  bifurcations.  This  investigation  illustrates  the  possible  rich 
dynamics  in  multibody  systems  with  closed  loop  structure  even  with  small 
number  of  degrees  of  freedom. 

*  This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative 
Program  under  grants  AFOSR-87-0073  and  AFOSR-90-0105  and  by  the  National  Science 
foundation’s  Engineering  Research  Centers  Program:  NSFD  CDR  8803012,  and  also  by  the 
Army  Research  Office  through  the  Mathematical  Sciences  Institute  of  Cornell  University. 
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1.  Introduction 

An  interesting  problem  in  the  mechanics  of  multibody  systems  is  to  find  relative 
equilibria.  Investigations  along  this  line  in  several  examples,  e.g.  planar  two-  or  three- 
body  system,  two-body  system  with  ball-in-socket  joints,  etc.,  have  been  reported  in  [2, 
6,9-11].  Increasing  sophistication  in  spacecraft  configurations  (e.g.  the  NASA  Flight 
Telerobotic  Servicer,  modular  space  station  designs,  etc.)  has  underscored  the  need  for 
a  better  understanding  of  the  dynamics  of  multibody  systems.  Systems  incorporating 
kinematic  loops,  as  does  the  four-bar  linkage,  are  particularly  challenging.  The  ideas  of  the 
present  paper  might  be  used  for  instance  in  bringing  a  modular  articulated  space  structure 
into  a  nonstandard  (e.g.  nonsymmetric)  equilibrium  shape  by  slowly  altering  a  bifurcation 
parameter,  say,  through  mass  shifting.  It  is  our  intention  in  this  paper  to  illustrate  the 
possibilities  for  such  a  control  strategy  by  focusing  on  the  bifurcation  phenomena  in  one 
example,  the  four-bar  linkage.  In  [12]  we  studied  the  geometry  and  dynamics  of  floating 
four-bar  linkages  which  may  have  a  general  structure.  Relevant  notations,  terminology, 
and  key  concepts  are  recalled  in  Appendix  A.  In  this  paper,  we  study  qualitatively  the 
relative  equilibria  of  a  floating  four-bar  linkage  with  symmetric  configurations. 

Referring  to  the  notations  in  Appendix  A,  a  floating  four-bar  linkage  is  of  symmetric 
type  if,  with  proper  labels  of  the  bars, 


mi  =  m  3 


doi  —  d03,  dio  =  d30,  di2  =  d32,  d2i  =  d23-  (1*1) 

In  other  words  it  can  form  symmetric  shapes  as  shown  in  Fig.  1.1.  It  is  not  hard  to  verify 
that  any  four-bar  linkage  which  satisfies  (1.1)  has  two  such  symmetric  configurations. 

In  this  paper  we  shall  show  how  the  relative  equilibria  of  this  system  -vary  with 
respect  to  certain  parameters:  the  positions  of  body-centers-of-mass.  We  shall  show  that 
the  two  symmetric  configurations  are  relative  equilibria  for  any  choice  of  the  parameters. 
By  varying  the  offset  of  the  center  of  mass  of  the  0-th  bar,  the  symmetric  relative  equilibria 
{ire  bifurcated  generically.  The  bifurcations  are  of  pitchfork  type.  For  different  choices  of 
position  of  body-center-of-mass  of  the  2-nd  bar,  this  pitchfork  bifurcation  can  be  either 
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Fig.  1.1  A  Symmetric  Configuration 

supercritical  or  subcritical.  To  avoid  tedious  calculations,  we  will  consider  a  simplified  case 
which  also  exhibits  these  bifurcation  phenomena.  In  fact,  the  verification  of  the  assertions 
in  this  paper  have  been  carried  out  both  numerically  and  symbolically.  For  symbolic 
analysis,  MACSYMA1  was  employed. 

As  in  [12],  the  theorem  of  Smale  has  been  used  to  determine  relative  equilibria  of 
the  system.  The  singularity  theory  [3,5]  has  been  applied  in  the  analysis  of  the  bifurcation 
problem. 

2.  Relative  equilibria  and  bifurcations 

We  recall  Smale’s  theorem  about  relative  equilibria  [8].  Consider  a  hamiltonian 
system  ( )  where  M  is  a  2n -dimensional  smooth  manifold,  w  is  a  symplectic 

1  MACSYMA  is  a  trademark  of  Symbolics  Inc.,  Cambridge,  Mass. 
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form  on  A/,  Xu  is  hamiltonian  vector  field  determined  by 


tXflW  =  dH 

for  some  smooth  hamiltonian  function  H  on  M .  Here,  we  consider  natural  mechanical 
systems,  so  M  =  T*Q ,  the  cotangent  bundle  of  the  configuration  space  Q,  of  the  system. 
Let  G  be  a  Lie  group  acting  on  Q  on  the  left, 

Q-.GxQ^Q 

( 9,q )  •-»  $,(«)  = 

and  hence  on  T*Q  by  cotangent  lift.  The  following  definition  of  a  relative  equilibrium  is 
standard  [1]. 

Definition  2.1:  Let  n.  be  the  flow  of  Xu  on  M .  Then  zt  6  M  is  a  relative  equilibrium 
if  FxH(ze)  is  a  stationary  motion,  i.e.  there  exists  £  €  Q  such  that 

FxBM  =  exp(t()(ze), 

where  Q  is  the  Lie  algebra  of  the  group  G . 

Remark:  A  physical  interpretation  of  relative  equilibria  is  that  if  the  dynamics  of  a  system 
is  rotationally  invariant,  the  dynamical  orbit  of  a  relative  equilibrium  appears  to  be  a  fixed 
point  for  an  observer  in  a  suitable  uniformly  rotating  coordinate  system. 

Given  a  simple  mechanical  system  with  symmetry,  ( Q,K,V,G ),  where  K  is  a  G- 
invariant  Riemannian  metric  on  Q,  V  :  Q  —+  R  is  a  G-invariant  potential  function,  we 
have  the  following  very  useful  theorem. 

Theorem  :  [8]  For  simple  mechanical  system  with  symmetry  (Q,K,  V,  G),  define 

V(  :  Q  -  R  :  ,  ~  V(t)  -  ±K((Q(qM<iM)  -  (2.1) 

for  each  (  6  Q ,  where  is  infinitesimal  generator  of  the  action  corresponding  to  £ .  Then 
zt  —  (qe,pe)  &  T*Q  is  a  relative  equilibrium  if  and  only  if  qt  is  a  critical  point  of  for 
some  £  €  G  and  pe  =  K ^((Q{qt))- 
Remark: 

(1).  The  function  Vj  is  called  augmented  potential  function. 
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(2).  It  can  be  shown  that,  for  a  given  £  £  £?,  has  the  symmetry, 


W,(*))  =  V((x) 

for  all  g  £  Gj  :=  {g  £  G|Ady£  =  £} .  If  G  is  abelian,  G ^  —  G.  If  the  action  <3?  is  free 
and  proper,  then  the  quotient  space,  Q/G ^ ,  is  a  smooth  manifold  and  ir^  :  Q  — ♦  Q/G ^  is 
a  submersion.  Thus  Vjt  induces  a  function  on  Q/G {  such  that 

=  Vj  0  7T{. 

(3) .  If  G  is  abelian,  it  can  be  shown  that  if  7r z(qe)  is  a  local  minimize!  of  , 
the  corresponding  zt  =  ( qe,pe )  is  a  stable  relative  equilibrium,  and  if  n^(qe)  is  a  local 
maximizer  of  ,  the  corresponding  zt  =  (qe,Pe)  is  an  unstable  relative  equilibrium.  We 
will  refer  to  Q/G $  as  the  shape  space ,  and  the  points  ^^(qe)  as  the  relative  equilibrium 

shapes. 

(4) .  One  can  also  use  the  amended  potential  in  the  sense  of  Smale  to  study  the 
relative  equilibria  (see  [7,8]).  | 

Now  we  return  to  the  problem  of  relative  equilibria  of  floating  four-bar  linkages  which 
admit  symmetric  configurations.  In  [12]  we  have  shown  that  a  floating  four-bar  linkage  is 
a  simple  mechanical  system  with  symmetry  with  group  S1 .  Moreover,  setting  Vr  to  be 
zero,  locally,  Vjt  is  of  the  form 

VKM  =  -erJe  (2.2) 

where  e  =  (1  1  1  I)71,  the  elements  of  the  4x4  matrix  J  are  functions  of  a  relative  angle, 
say,  #10  =  #1  ~~  $0  •  Although  it  is  not  easy  to  find  analytically  the  critical  points  of  (2.2), 
for  a  particular  example  one  can  easily  do  this  numerically  since  Vjc  is  only  a  one  variable 
function.  Unlike  the  planar  two  body  problem  [9]  for  which  the  dimension  of  the  shape 
space  is  also  one,  here  the  function  VJc  has  many  parameters,  even  under  the  conditions 
of  (1.1).  A  natural  question  is  to  determine  how  these  parameters  affect  the  relative 
equilibria,  e.g.  their  numbers  and  location  on  shape  space,  etc..  Of  course,  it  is  difficult  to 
answer  this  question  for  completely  arbitrary  choice  of  parameters.  However,  by  leaving  one 
particular  parameter  free  such  that  the  assembly  preserves  its  symmetric  configurations  and 
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freezing  all  other  parameters,  one  still  can  observe  a  nontrivial  bifurcation  phenomenon. 
To  illustrate  this  we  consider  an  example. 

Example:  Let  us  choose  the  parameters  as  follows. 

m0  =  1,  mi  =  1,  m2  =  1,  m3  =  1; 

do3=(— 2,Ao),  doi  =  (— 2, A0),  -  dio  =  (— 1.5,  — 1),  di2  =  (1.5,  — 1), 

d21  =  (-1.5,  -1.4),  d23  =  (1.5,  -1.4),  d32  =  (-1.5,  -1),  d30  =  (1.5,  -1). 

Now  the  assembly  has  non-Grashof  structure.  We  know  that  6io  cannot  be  used  as  a 
global  coordinate  system  since  for  each  tfio  there  are  two  points  in  configuration  space, 
which  has  been  shown  to  be  a  torus  [4].  One  point  corresponds  to  leading  form.  Another 
one  corresponds  to  lagging  form  (see  Appendix  A).  Nevertheless,  one  can  see  that  6\q 
taken  together  with  the  sign  of  sm(023)  can  separate  these  forms.  For  this  reason,  we 
define  a  new  variable  as  follows.  Assume  the  0-th  bar  is  the  longest  bar.  Let  a  >  0  be 
the  maximum  attainable  interior  angle  between  0-th  bar  and  1-st  bar.  Then  6  is  defined 
a s  follows: 

—  y,  if  sm(023)  >  0  (leading  form); 

*19-*-.?.  if  sin(023 )  <  0  (lagging  form). 

It  is  clear  that  if  we  identify  —  n  and  n  to  be  the  same  point  on  S1 ,  6  parameterizes  51 . 

Now,  using  Vj,  for  any  Ao  one  can  find  relative  equilibria  (#io)t  for  both  leading 
form  and  lagging  form,  and  hence  one  can  find  8.  As  Ao  changes  from  —00  to  +00, 
one  can  plot  a  diagram  with  A0  and  8  as  the  parameter.  Fig.  2.1  shows  the  result,  in 
which  solid  dots  represent  stable  relative  equilibria,  small  circles  represent  unstable  relative 
equilibria. 

From  this  example  one  can  make  the  following  empirical  observations: 

(1) .  There  are  two  unbounded  symmetric  branches  on  the  diagram  and  and 
these  branches  are  bifurcated  at  some  points.  The  bifurcations  appear  to  be  pitchfork 
bifurcations. 

(2) .  Almost  any  value  of  6  can  be  relative  equilibrium  for  a  particular  A0.  In  other 
words  the  bifurcation  diagram  is  connected  globally. 
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X0 

Fig.  2.1  Bifurcation  diagram:  an  example 
(3).  The  number  of  relative  equilibria  can  be  two,  six  and  ten. 

The  first  observation,  which  relates  to  the  local  bifurcation  problem,  is  what  we  will 
concentrate  on  in  the  rest  of  this  paper.  The  others  will  be  discussed  in  section  4. 

3.  The  recognition  problem  in  bifurcation  theory 

Our  problem  raised  in  the  previous  section  can  be  reduced  to  showing  how  the 
solutions  x  of  a  single  scalar  equation 

9{xi  A)  =  0  (3.1) 

change  with  the  parameter  A;  or,  more  precisely,  what  type  of  bifurcation  occurs  with 
parameter  A.  Without  loss  of  generality,  one  can  assume  <7(0,0)  =  0.  Moreover,  we 
assume  g  :  R  X  R  — *  R  is  smooth.  This  is  one  of  the  standard  local  (static)  bifurcation 
problems  with  one  state  variable,  called  the  recognition  problem.  As  with  many  bifurcation 
problems,  this  problem  can  be  solved  successfully  through  singularity  theory,  in  which  the 
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related  issue  is  called  finite  determinacy.  Let  z  —  (2,  A).  Near  origin,  the  function  g  can 
be  written  as 

»w=  E  ^(^)“9(0)z”+  E  M*)*".  (3.2) 

H<fc+i  '  *  |o|=t+i 

for  some  smooth  functions  aa  defined  in  a  neighborhood  of  the  origin.  Here  we  used  the 
conventions  with  multi-indices: 


|q|  =  au  +  at2,  a!  =  (qi)!(q2)!, 


.<*1 


A  key  question  is  what  terms  in  Eq.(3.2)  can  be  ignored  such  that  the  values  of  coefficients 
of  remaining  terms  can  be  used  to  determine  the  qualitative  behavior  of  the  original 
equation  (3.1),  for  example,  the  variation  in  number  of  solutions.  Singularity  theory  solves 
this  problem  by  finding  a  suitable  change  of  coordinates  such  that  function  g  is  equivalent 
to  a  standard  model  h ,  called  normal  form.  A  precise  definition  is  given  as  follows  (see 
[5])- 


Definition  3.1:  Two  smooth  mappings  g,h  :  R  x  R  — ►  R  defined  near  the  origin  are 
equivalent  if  there  exist  a  local  diffeomorphism  of  R2 ,  (r,A)  (A'’(r,  A),  A(A))  at  the 

origin  and  a  nonzero  function  S(x,  A),  such  that 


g(x,  A)  =  S(r,  \)h(X(x,  A),A(A))  (3.3) 

where  Xx(0,0)  >  0  and  A'(0)  >0.  If  A  =  A,  g  and  h  are  strongly  equivalent. 

From  this  definition  we  see  that,  since  S(x,  A)  is  nonzero,  the  solution  of  g(x,  A)  =  0 
and  h(X,  A)  =  0  are  the  same  in  the  sense  of  diffeomorphism.  From  this  point  of  view,  by 
means  of  singularity  theory  one  can  show  why  and  what  the  high-order  terms  in  Eq.(3.2) 
do  not  effect  the  qualitative  behavior  of  equation  g(x,  A)  =  0.  It  should  be  noticed  that 
although  this  method  does  not  tell  us  how  to  derive  an  appropriate  normal  form  h ,  for 
most  physical  problems,  such  as  the  one  considered  in  this  paper,  it  is  not  hard  to  pick  up 
some  of  the  candidates  from  a  large  number  of  known  simple  polynomials  of  x  and  A,  or 
the  model  of  normal  forms  which  have  standard  bifurcation  diagrams.  This  is  in  essence 
the  spirit  of  application  of  singularity  theory  to  a  physical  problem. 
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Without  considering  detailed  issues  of  singularity  theory  which  are  applicable  to 
bifurcation  problems,  we  directly  give  the  following  result  which  will  be  used  in  the  next 
section.  For  details  see  [5],  chapter  2.  First  we  need  the  concept  of  germs. 

Definition  3.2:  Two  smooth  functions  defined  near  the  origin  are  equivalent  as  germs  if 
there  is  some  neighborhood  of  the  origin  on  which  they  coincide.  Let  Sz  \  denote  the  set 
of  equivalence  classes  of  such  functions.  The  elements  in  Sx\  are  called  germs. 

Lemma:  A  germ  g  £  £Zt\  is  strongly  equivalent  to 

exk  4-  6Ar  (3-4) 

for  k  >  2  if  and  only  if  at  x  —  A  =  0 

s  =  A9  =  ...  =  <£)*-’»  =  |^=0  (3'5)“ 


smd 


«  =  S9n(g^)k9 


c  d  d 


(3.5)6 


Remark:  When  k  =  3  the  normal  form  (3.4)  provides  a  pitchfork  bifurcation.  From  this 
lemma,  so  is  g  if  (3.5)  holds.  It  is  easy  to  show  that  if  e6  >  0,  the  pitchfork  bifurcation  is 
subcritical;  if  eS  <  0,  the  pitchfork  bifurcation  is  supercritical. 

4.  Bifurcations  at  symmetric  configuration 

As  we  have  seen,  the  function  of  a  four-bar  linkage  is  a  multiple  parameter 
function.  One  might  expect  very  complicated  bifurcation  features  with  respect  to  these 
parameters.  Here,  instead  of  considering  a  general  structure,  we  study  a  special  assembly 
which  has  symmetric  configuration  as  defined  in  section  2.  To  avoid  too  many  tedious 
calculations  we  particularly  choose  the  parameters  of  the  assembly  as  follows. 


m0  =  m  1  =  m2  =  m3  =  1; 

do3  =  (— d0,A0),  d0i  =  (do,Ao),  dio  =  (—1,0),  dj2  =  (1,0)  (4.1) 

d2i  =  (— <^2,  A2),  d23  =  (d2,  A2),  d32  =  (  — 1,0),  d3o  =  (l,0). 
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where  d0  and  d2  are  fixed  and  d0  >  d2  >  0,  A0,A2  £  R-  Moreover,  we  consider  the 
non-Grashof  case  only,  i.e. 

do  +  d2  +  2  —  2(max{do, c?2, 1}  +  mm{d0,  d2, 1})  <  0.  (4.2) 

Figure  4.1  shows  two  symmetric  configurations  for  above  choice  of  parameters.  We 
will  see  that  although  only  two  parameters  Ao  and  A2  are  left  to  be  free,  the  bifurcation 
features  with  respect  to  these  parameters  are  still  informative. 


I  mo 


Fig.  4.1  Symmetric  Configurations 

The  function  \\  now  has  the  following  form: 

=  i(2A2sm(032)  +  d2cos{d32)  ~  cos(83i)  +  d0cos(63o) 
+2X2sin(62i)  +  d2cos(821 )  -  d0d2cos(d20)  +  docos(0io)) 
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(4.3) 


-f^j-(2sm(03o)  +  X2cos(020)  -  2<sm(0lo))  +  C 

where  C  is  a  constant  determined  by  d;2,  m*  and  the  moments  of  inertia  of  the  bodies, 
Ii.  For  i,j  =  0,  1,  2,  3,  6,j  =  9%  —  8j  and  6i  satisfy  the  constraint  equations 


d0  +  cos(0lo)  +  d2cos(620)  +  cos(630)  =  0 


(4.4)a 


sin{9 io )  +  d2sin(920)  +  sin(^0)  =  0. 


(4.4)6 


In  the  following,  at  symmetric  configuration,  the  variables  will  be  denoted  by 
superscript  “s”  (say,  #i0),  the  formulas  will  be  denoted  by  (say,  /(#io)|»)-  As 

shown  in  the  example  in  section  2,  the  bifurcation  diagram  of  relative  equilibria  will  be 
parameterized  by  (6,  Ao). 

Theorem  4.1:  For  a  floating  four-bar  linkage  with  parameters  shown  in  Eq.(4-l),  the 
bifurcation  diagram  of  relative  equilibria  has  the  following  properties: 

(1) .  There  are  two  infinite  branches  in  the  diagram,  one  corresponds  to  8S10  in  the 
leading  form,  another  one  corresponds  to  in  the  lagging  form.  We  refer  to  these  as  the 
symmetric  branches. 

(2) .  There  exists  a  constant  A2  such  that  no  bifurcation  occurs  on  the  symmetric 
branch  of  leading  form  if  X2  =  A2 ;  and,  no  bifurcation  occurs  on  the  symmetric  branch  of 
lagging  form  if  \2  =  —  A2 . 

(S).  On  the  symmetric  branch  of  leading  (lagging)  form,  if  A2  <  A2  (  — A2),  there 
exists  a  constant  C\  (03)  such  that  the  relative  equilibria  are  stable  for  Ao  <  C\  (03),  unstable 
for  A0  >  ci  (C3),  bifurcated  for  A0  =  Ci  (C3);  one  the  other  hand,  if  \2  >  A2  (  — A2),  there 
exists  a  constant  c2  (C4)  such  that  the  relative  equilibria  are  unstable  for  A0  <  c2  (C4), 
stable  for  Ao  >  c2  (c4),  bifurcated  for  Ao  =  C\  (c4). 

(4).  Assume  X2  ^  ±A2  •  Let 


=  sgn 


ieiAj  4-  C2 A2 i^3 
€4  A2±e5 


(4.5) 


and 


6*  =  sgn(6i  A2±<52) 


(4.6) 
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where  e ,  and  Si  are  constants  which  are  determined  by  do  and  di .  Then,  the  bifurcation 
on  the  symmetric  branch  of  leading  form  will  be  supercritical  pitchfork  if  <  0  .It  will 

be  subcritical  pitchfork  if  e+6+  >  0.  Similarly,  the  bifurcation  on  the  symmetric  branch  of 
lagging  form  will  be  supercritical  pitchfork  if  e~6~  <0.  It  will  be  subcritical  pitchfork  if 
e~6~  >  0. 

Remark: 

(a) .  Based  on  the  techniques  of  bifurcation  theory  shown  in  section  3,  the  proof  of  the 
above  assertions  is  very  elementary.  However,  two  aspects  have  to  be  considered  before  the 
bifurcation  theory  is  applied.  First,  one  has  to  determine  the  existence  of  the  bifurcation 
and  where  the  branches  are  bifurcated.  Second,  as  one  can  see,  it  is  not  easy  to  write 
down  explicitly  the  function  as  a  function  of  one  variable.  Therefore,  when  applying 
the  techniques  of  bifurcation  theory  which  involves  as  high  as  fourth  order  derivatives  of 
the  function  Vj ,  one  should  consider  the  constraint  equations  simultaneously. 

(b) .  Although  the  proof  of  the  above  assertions  is  elementary,  it  requires  a  large 
effort  in  calculations.  We  used  MACSYMA  to  handle  these  computations  (see  listings  at 
the  end  of  the  paper).  In  the  following,  we  only  give  a  sketch  of  the  proof. 

Proof  of  theorem  4.1:  Note  that  the  function  V, j  can  be  written  as  a  function  of  relative 
angles  #10,  #20,  and  #30  ,  which  are  related  through  constraint  equations  (4.4).  From.  (4.4), 
we  can  consider  (locally)  #20  and  $30  as  the  functions  of  0io-  Again,  from  (4.4)  one  can 
generate  the  quantities  -q^\»  and  |,  for  any  positive  integer  i.  Moreover,  from  Fig. 
4.1  it  is  easy  to  see  that  at  symmetric  configuration 


620  —  ft 


and 


*;•>  = 


- 0 s 


30- 


(4.7) 


With  above  considerations,  one  can  have  closed  form  expressions  for 


For  instance, 


80 


22. 


80 10 

802d  I  —  1 

801O  I*  ~ 


=  j;cos{6{  0), 


8'0 


22. 


and 


8'0 


2SL 


<3*lo 


(4.S) 


and 


2co$(0‘io) 

d7»in(0[o) 


(d2  -  2cos(9’0)) 


(4.9) 
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To  prove  assertion  (1)  in  the  statement  of  the  theorem,  one  needs  to  show  that  at 
symmetric  configuration,  the  equation 


dVj 

dd\o 


=  0 


does  not  depend  on  Ao .  Concentrating  on  the  term  involving  A0  in  and  applying 

Eq.(4.7)  and  (4.8),  one  can  show  that  the  first  derivative  of  that  term  with  respect  to 

dv 

0IO  at  symmetric  configuration  is  zero.  Since  the  rest  of  the  terms  of  =  0  are  still 

functions  of  ,  one  can  see  two  infinite  symmetric  branches  in  the  bifurcation  diagram 
for  two  different  0^o .  Assertion  (1)  is  thus  proved. 

Applying  (4.8)  and  (4.9),  one  can  show  that  the  second  derivative  of  the  function  \\ 
at  symmetric  configuration  has  the  form 


s2% ,  _  yn  _  ,<h-2 <W(g;„)  <W(a;0) 

ae2J'  ~  ae'i0]‘  M{  d2sm(e>u)  2  4 


(4.10) 


where  II  is  the  summation  of  the  terms  not  involving  Aq  in  .  It  is  obvious  that  when 


x  _  y  A  2 d2cos2(0’iO)  -  (P2 
'2  cos2(0’o)sin(0* 0)’ 


(4.11) 


^2  y  '  ^2  y 

will  not  depend  on  A0 .  One  can  also  show  that  with  Eq.(4.11),  7^  0 

under  assumptions  (4.1)  and  (4.2).  This  means  that  bifurcation  may  not  occur  on  either 
symmetric  branch  of  leading  form  or  symmetric  branch  of  lagging  form.  Note  that  on 
these  different  forms  cos(0[o )  has  the  same  value,  sin(0{o )  has  the  same  absolute  value 
but  different  sign.  (See  Fig.  4.1)  Thus,  assertion  (2)  is  proved. 

As  we  have  known  earlier,  the  stability  of  relative  equilibria  depends  on  the  sign  of 

^2 y  q2 y 

'a^r-U-  From  Eq.(4.10)  we  see  that  is  a  linear  function  of  Ao .  Using  the  \2  in 

(4.11),  the  proof  of  (3)  is  straight  forward. 

To  prove  assertion  (4),  we  apply  the  bifurcation  theory  mentioned  in  section  (3).  Let 
AJ  denote  ct  in  assertion  (3)  for  some  suitable  i.  One  can  show  that  at  (ffio,  AjJ), 


Moreover 


(4.13) 


1*  +  ^2 -^2 i^3 

a^0(l0’  o)~  64A2±e5 

and 

Snj^Wc  A5)  =  Mj±«2  (4.14) 

for  some  constants  e,  and  <5, ,  where  “  4-  ”  corresponds  to  leading  form,  “  —  ”  corresponds 
to  lagging  form.  See  Appendix  B  for  expressions  of  and  Si.  Since  (4.13)  and  (4.14) 
are  not  zero  in  general,  applying  the  lemma  in  section  3,  we  can  say  is  strongly 

equivalent  to  the  normal  form  of  pitchfork  bifurcation.  In  addition,  the  type  of  pitchfork 
bifurcation  depends  on  the  sign  of  (4.13)  and  (4.14).  The  assertion  (4)  is  proved. 

1 

Remark:  (i).  The  condition  of  A2  AJ  guarantees  that  (4.14)  and  the  denominator  of 
(4.13)  are  not  zero.  (ii).  In  general  the  bifurcation  changes  from  a  supercritical  one  to 
subcritical  one  at  the  roots  of  numerator  of  (4.13). 

Example:  To  see  how  A2  changes  the  bifurcation  diagram  with  parameter  A0  we  give 
following  example.  We  will  concentrate  on  the  symmetric  branch  with  respect  to  leading 
form.  Let  d0  —  2  and  d2  =  1 .  Then  e  and  6  have  the  following  form 

,  \\  +  6.938A2  -  3.623 

y  1.873  -  0.051A2 

6+  =50n(l.O25-O.O28A2) 


Then 


+6+  _  I  -1,  if  -7.426  <  A2  <  0.488; 
~  \+l,  otherwise. 


Note  that  the  region  for  A2  is  an  approximation.  So  we  can  say  that,  when  A2  € 
(—7.426,  0.488),  the  pitchfork  bifurcation  is  supercritical.  Otherwise,  it  is  subcritical. 
Fig  4.2  shows  this  result. 


Before  closing  this  section,  we  would  like  to  make  some  further  remarks  regarding 
this  paper. 

(i).  Although  our  discussion  only  concentrated  on  a  structure  of  non-Grashof  type, 
a  version  of  theorem  4.1  also  holds  for  the  Grashof  case. 
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(ii) .  Up  to  now  we  have  understood  the  phenomenon  of  bifurcation  on  the  symmetric 
branches.  The  globed  analysis  of  the  bifurcations  is  in  progress.  This  involves  massive 
symbolic  computations.  However,  as  shown  in  the  example  in  section  2  ancl  other 
simulations,  one  can  numerically  determine  a  global  bifurcation  diagram.  A  large  body 
of  such  numerical  simulations  show  that  the  branches  in  the  bifurcation  diagram  are 
connected.  In  other  words,  for  any  point  in  shape  space,  there  is  a  finite  Ao  for  which 
that  shape  determines  a  relative  equilibrium.  This  is  consistent  with  the  linearity  of  VJt  in 
Aq  .  This  property  provides  a  possibility  to  control  the  attitude  of  a  space  structure  with  a 
closed  kinematic  chain  by  simply  changing  the  position  of  the  center  of  mass  of  some  bars. 

(iii) .  Our  results  in  this  paper  rely  on  some  ideal  conditions,  for  instance,  the 
symmetry  conditions  (1.1),  and  the  absence  of  external  and  internal  disturbances.  One 
may  ask  what  will  happen  when  these  conditions  are  violated.  The  answer  to  this  question 
will  be  related  to  the  notion  of  universal  unfolding  in  bifurcation  theory.  We  will  explore 
this  in  a  further  paper.  Numerical  results  show  that  it  is  possible  to  use  the  unfolding 
property  to  control  the  shape  of  the  structure  near  the  bifurcation  point. 

5.  Conclusion 

We  have  investigated  the  local  and  global  structure  of  bifurcations  of  relative 
equilibria  in  a  model  multibody  problem  -  the  floating  four-bar  linkage.  The  observed 
global  connectivity  of  the  bifurcation  diagram  suggests  that  it  might  be  possible  to 
change  the  attitude/orientation/shape  of  such  a  system  by  merely  tuning  an  appropriate 
parameter,  e.g.  a  mass  offset.  The  precise  dynamics  associated  to  such  a  control  strategy 
would,  of  course,  depend  on  the  role  of  internal  dissipation  and  consequent  change  of  a 
stable  center  into  a  focus.  Further  work  on  the  phase  portrait  and  domains  of  attraction 
is  under  way  and  we  hope  to  report  on  these  matters  in  a  future  paper. 

Since  we  employed  the  symbol  manipulation  language  MACSYMA  for  our  local 
analysis,  we  have  included  a  listing  of  the  relevant  MACSYMA  program  at  the  end  of 
this  paper,  to  assist  the  reviewer. 
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Appendix  A 

A.l.  Notations  and  Geometric  Constraints 

The  structure  of  a  closed  floating  four-bar  linkage  is  represented  in  Fig.  A.l.  The 
bars  are  labeled  clockwise  from  0  to  3  as  shown.  We  define  the  following  quantities. 


V 


Fig.  A.l  The  general  structure  of  four-bar  linkage 


d  ij 


R(0,) 


the  vector  of  hinge  point  which  connects  i-th  bar  with  j -th  bar 
relative  to  a  body-fixed  frame  with  origin  at  the  center  of  mass 
of  the  i-th  body; 


the  rotation  through  angle  0,  giving  the  orientation  of  i-th  bar 


relative  to  the  inertial  space; 


-  m)  = 


(  cos(6t ) 
ysin(6i) 


cos(0,)  J  ’ 


joint  rotation  between  i  -th  and  j  -th  bar, 


R(6ij)  =  R{G,  -  Oj)  = 


IS 


the  distance  between  the  joints  on  t-th  bar,  (or  ’’length”  of  t-th 
bar);  all  /,  >  0;  /,  =  ||dl)t+1  -  dif<_i||; 


U 


With  above  notations,  the  loop  constraint  can  be  represented  as 

3 

]T.R(0,)(d,it+1  -  d^i-j)  =  0, 
i= 0 

where  we  adopt  the  convention  that  dl4  =  d,to  and  d4);  =  do,,  . 

A. 2.  The  Configuration  Space 
Let 

s  =  length  of  shortest  bar 
l  =  length  of  longest  bar 
p,q  =  lengths  of  intermediate  bars. 

Assume  that  the  observer  is  at  the  center  of  mass  of  the  system.  Then  we  have  following 
results.  (See  [12]) 

If 

s  +  1 <  p+  q 

the  configuration  space  is  isomorphic  to  two  separated  tori.  Under  this  condition  the 
assembly  is  called  Grashof  type.  If 


s  -fi  l  >  p  -fi  q 

the  configuration  space  is  isomorphic  to  one  torus.  Under  this  condition  the  assembly  is 
called  non- Grashof  type.  If 

s  +  l  =  p  +  q 

the  configuration  space  is  not  a  smooth  manifold. 

Let  /j  =  s  and  define  the  orientation  of  bars  as  shown  in  Fig  A.l.  Then  we  say  that 
the  structure  is  in  leading  form  if  sin{8 3  —  $2)  <  0,  in  lagging  form  if  atn(^  —  02 )  >  0. 
From  this  definition,  we  see  that  for  the  Grashof  case,  each  component  of  the  configuration 
space  corresponds  to  a  distinct  form.  Two  such  forms  cannot  be  deformed  into  each  other. 
But  for  the  non-Grashof  case,  they  can  be  continuously  deformed  on  the  torus. 
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Appendix  B 

The  expression  of  and  6{  in  the  theorem  4.1  are  of  the  following  form. 

C\  =  — 12do(d2  —  do)2(dl  —  dod%  —  dod2  ~  2c?2  +  d0  —  4do), 
e2  =  48 d20d2((d2  -  d0)2  -  1  )(dl  -  d2)y/4-(d0-d2)2; 
e3  =  4&dldl(d2  —  do)2(dl  —  dodl  —  d\d2  —  4  d2  +  dl  —  2do); 
c4  =  ^((^2  —  do)2  —  4)(d2  —  do)2  y/ 4  —  (do  —  d2)2; 
e5  =  2 d*{{d2  -  do)2  -  4){(d2  -  do)3  -  4 d2). 

Let 

A  =  4d]  y/4  -  (d0  -  d2)2. 

Then  _ _ 

6\  =  —  ~^d2  —  do)2  y/ 4  —  (do  —  d2)'1\ 

62  —  —  ^d2{d\  +  Zdod\  —  3do<f2  +  4d2  +  do). 


20 


REFERENCES 


[1]  Abraham,  R.  &  J.  E.  Marsden,  Foundations  of  Mechanics  (Second  edition),  Revised, 
Enlarged,  Reset,  Benjamin/Cummings,  Reading,  1978. 

[2]  Baillieul,  J.,  “An  enumerative  theory  of  equilibrium  rotations  for  planar  kinematic 
chains,”  in  Dynamics  and  Control  of  Multibody  Systems,  J.  E.  Marsden,  P.  S. 
Krishnaprasad  &:  J.  C.  Simo,  ed.,  AMS  Contemporary  Mathematics  #Vol.  97,  1988, 
1-10. 

[3]  Brocker,  T.,  Differentiable  Germs  and  Catastrophes,  Cambridge  University  Press,  1975. 

[4]  Gibson,  C.  G.  &  P.  E.  Newstead,  “On  the  Geometry  of  the  Planar  4-Bar  Mechanism,” 
Acta  Applicandae  Math.  7  (1986),  113-135. 

[5]  Golubitsky,  M.  &  D.  G.  Schaeffer,  Singularities  and  Groups  in  Bifurcation  Theory, 
Volume  I,  Springer- Verlag,  1985. 

[6]  Patrick,  G.  W.,  “The  dynamics  of  two  coupled  rigid  bodies  in  three  space,”  in  Dynamics 
and  Control  of  Multibody  Systems,  J.  E.  Marsden,  P.  S.  Krishnaprasad  J.  C.  Simo, 
ed.,  AMS  Contemporary  Mathematics  #Vol.  97,  1988,  315-336. 

[7]  Simo,  J.  C.,  T.  Posbergh  &  J.  E.  Marsden,  Stability  of  Coupled  Rigid  Body  and 
Geometrically  Exact  Rods:  Block  Diagonalization  and  the  Energy-Momentum  Method , 
Stanford  University,  Division  of  Applied  Mechanics  Technical  Report  89-2,  1990,  (to 
appear  in  Physics  Report). 

[8]  Smale,  S.,  “Topology  and  Mechanics,  I,  II,”  Invent.  Math.  11  (1970),  45-64. 

[9]  Sreenath,  N.,  “Modeling  and  Control  of  Multibody  Systems,”  University  of  Maryland, 
College  Park,  Ph.D.  Thesis,  1987,  also,  Systems  Research  Center  Technical  Report  SRC 
TR87-163. 

[10]  Sreenath,  N.,  Y.  G.  Oh,  P.  S.  Krishnaprasad  &  J.  E.  M.  (1988),  “The  Dynamics  of 
Coupled  Planar  Rigid  Bodies  Part  I:  Reduction,  Equilibria  &:  Stability,”  in  Dynamics 
&:  Stability  of  Systems  3  #1&2,  25-49. 


21 


[n]  WANG,  L.  -S.  &  P.  S.  KRISHNAPRASAD,  “Relative  Equilibria  of  Two  Rigid  Bodies 
connected  by  a  Ball-in-Socket  Joint,”  Proc.  of  the  1080  IEEE  Conference  on  Decision 
and  Control, Y ol.  1  (Dec.  1989),  692-697. 

[12]  Yang,  R.  &  P.  S.  Krislmaprasad,  “On  the  Dynamics  of  Floating  Four-bar  Linkages,” 
Proceedings  of  IEEE  Conference  on  Decision  and  Control,  V ol.  2  (Dec.  1989.),  1632- 
1635. 


00 


The  MACSYMA  Program 
to  Determine  Local  Bifurcations 


/*******************************************/ 

/*  BIFURCATIONS  OF  RELATIVE  EQUILIBRIA  */ 
/*******************************************/ 

apply (writefile, [bifur_out] ) ; 

depends ( [th_20,  th_30],  th_10) ; 

/*********************************************/ 
/*  CONSTRAINT  EQUATIONS  */ 

/*********************************************/ 
f 1 :  d_0+cos (th_10>  +d_2*cos (th_20) +cos (th_30) $ 

£2 :  ain (th_10) +d_2*sin (th_20) +sin (th_30) $ 

^*****«***********************  ****************/ 
/*  THE  FIRST  DERIVATIVE  OF  CONSTRAINT  */ 

/*  EQUATIONS  W.R.T.  th_10  */ 

/*********************************************/ 
dlflrdiff (fl,  th_10,  1) $ 
dl£2 :dif f (f 2,  th_10,  1) $ 

/*********************************************/ 


/*  THE  EVALUATION  OF  */ 
/*  THE  FIRST  DERIVATIVE  OF  CONSTRAINT  */ 
/*  EQUATIONS  W.R.T.  th_10  AT  */ 
/*  SYMMETRIC  CONFIGURATIONS  */ 


/********«******»*****************************/ 
dlfls : ev (dlf 1 ,  diff (th_20,  th_10,  l)«dls_2, 
dif £ (th_30,  th_10,  l)-dls_3, 
th_20-%pi,  th_30— ths_10,  th_10-ths_10) $ 
dlf2s:ev  (dlf2,  diff  (th_20,  th_10,  l)»*dls_2, 
diff (th_30,  th_10,  l)-dls_3, 
th_20— %pi,  th_30— ths_10f  th_10“ths_10) $ 

/*********************************************/ 

/*  SOLVE  Dt h_2 0 / Dt h_l 0  AND  Dth_30/Dth_10  */ 

/*********************************************/ 
solve ( [dlf Is,  dlf2s],  [dls_2,  dls_3] ),  globalsolve : true ; 

/*********************************************/ 

/*  THE  SECOND  DERIVATIVE  OF  CONSTRAINT  */ 

/*  EQUATIONS  W.R.T.  th_10  */ 

/*********************************************/ 
d2f 1 : dif f (dlf 1 ,  th_10,  1) $ 
d2  f  2 : di f  f  (dlf2,  th_10,  1)$ 

/a********************************************/ 


/*  THE  EVALUATION  OF  */ 
/*  THE  SECOND  DERIVATIVE  OF  CONSTRAINT  */ 
/*  EQUATIONS  W.R.T.  th_10  AT  */ 
/*  SYMMETRIC  CONFIGURATIONS  */ 


/*  *******★**★*★*★★•**★★★★*★★**•*★★*★*★*★*■*★*★*★★/ 

d2fla :ev (d2f 1,  diff(th_20,  th_10,  l)«dls_2,  diff(th_30,  th_10,  l)-dls_3, 
diff (th_20,  th_10,  2)-d2s_2,  diff(th_30,  th_10,  2)-d2s_3, 
th_20-%pi,  th_30— ths_10,  th_10-ths_10) $ 
d2f 2s : ev  (d2f 2,  diff(th_20,  th_10,  l)-dls_2,  diff(th_30,  th_10,  l)-dls_3, 
diff (th_20,  th_10,  2)-d2s_2,  diff(th_30,  th_10,  2)-d2o_3, 
th_20— %pi,  th_30—ths_10,  th_10**ths__10)  $ 

/*************************************************/ 

/*  SOLVE  DDth_20/DDth_10  AND  DDth_30/DDth_10  */ 
/****«********************************************/ 
solve ( [d2f Is,  d2f2s],  [d2s_2,  d2s_3 ]), globalsolve : true ; 

d2s_2 : factor (d2 s_2 ) ; 
d2s  3 : factor (d2s_3 )  ; 


/  ★  ik  *  *  *  1  &***★★*★★***★*****★★**★★★*★★★*★★★★****★★/ 

/*  THE  THIRD  DERIVATIVE  OF  CONSTRAINT  */ 

/*  EQUATIONS  N.R.T.  th_10  */ 

^*********************************************/ 
d3fl :dif f (d2f 1,  th_10,  1)$ 
d3f2:d:Lff  (d2f2,  th_10,  1)$ 


/* 

THE  EVALUATION  OF 

*/ 

/* 

THE  THIRD 

DERIVATIVE  OF  CONSTRAINT 

*/ 

/* 

EQUATIONS 

W.R.T.  th_10  AT 

*/ 

/* 

SYMMETRIC 

CONFIGURATIONS 

*/ 

/******  *★★*****★*★*★★*★★*★**★★★***★★★★★*★★★♦*★/ 
d3fla :«v (d3£l,  diff(th_20,  th_10,  l)«dls_2,  diff(th_30,  th_10, 
dif f (th_20,  th_10,  2)-d2s_2,  diff(th_30,  th_10, 
diff (th_20,  th_10,  3)-d3s_2,  diff(th_30,  th_10, 
th_20-%pi,  th_30— ths_10,  th_10— ths_10) $ 
d3£2s :ev (d3f2,  diff(th_20,  th_10,  l)-dls_2,  diff(th_30,  th_10, 
diff (th_20,  th_10,  2)-d2s_2,  diff(th_30,  th_10, 
diff (th_20,  th_10,  3)-d3s_2,  diff(th_30,  th_10, 
th_20“%pi/  th_30“-ths_10,  th_10«“ths_10)  $ 


1)  «*dls_3, 

2) -d2s_3, 

3)  “d3s_3, 

1) -dls_J3, 

2) -d2s_3, 

3) -d3s  3, 


/♦it***************************************************/ 

/*  SOLVE  DDDt h_2 0 /DDDt h_l 0  AND  DDDth_30/DDDth_10  */ 
/*****************************************************/ 
solve ( [d3fls,  d3f2s],  [d3s_2,  d3s_3] ) , globalsolveitrue; 
d3s_2 : factor (d3s_2) ; 
d3s  3: factor (d3s  3); 


/*********************************************/ 
/*  THE  FOURTH  DERIVATIVE  OF  CONSTRAINT  */ 

/*  EQUATIONS  W.R.T.  th_10  */ 

/*********************************************/ 
d4fl :diff (d3fl,  th_10,  1)$ 
d4£2 : dif f (d3f2,  th_10,  1)$ 


/* 

THE  EVALUATION  OF 

*/ 

/* 

THE  FOURTH  DERIVATIVE 

OF  CONSTRAINT 

*/ 

/* 

EQUATIONS  W.R.T.  th_10 

AT 

*/ 

/* 

SYMMETRIC  CONFIGURATIONS 

*/ 

/*********************************************/ 
d4fls:ev(d4fl,  diff(th_20,  th_10,  l)-dls_2,  diff(th_30, 

diff (th_20,  th_10,  2)-d2s_2,  dif f (th_30, 

diff (th_20,  th_10,  3)-d3s_2,  diff(th_30, 

diff (th_20,  th_10,  4)-d4s_2,  diff(th_30, 

th_20«%pi,  th_30— ths_10,  th_10«ths_10) $ 
d4£2s :ev (d4f2,  diff(th_20,  th_10,  l)-dls_2,  diff (th_30, 

diff (th_20,  th_10,  2)-d2s_2,  diff(th_30, 

diff (th_20,  th_10,  3)-d3s_2f  diff(th_30, 

diff (th_20,  t h_l 0 ,  4)-d4s_2,  diff(th_30, 

th_20-%pi,  th_30— ths_10,  th_10**ths_10) $ 


th_10, 
t  h_l 0 , 
th_10, 
th_10, 

th_10, 
th_10, 
th_10, 
th  10, 


1) -dls_3, 

2) «d2s_3, 

3)  -d3s_3, 

4) -d4s_3, 

1) -dls_3, 

2) -d2s_3, 

3) -d3s_3, 

4)  -d4s_3. 


/■ink****  *********  +  *■********************  +  ****************  +  **/ 

/*  SOLVE  DDDDt h_2 0 /DDDDt h_l 0  AND  DDDDt h_3 0 / DDDDt h_l 0  */ 

/*«***  **+*++*+********  +  ***********************************/ 
solve ( [d4f Is,  d4f 2s] ,  [d4s_2,  d4s_3] )  ,  globalsolve :true; 

d4o_2 : factor (d4s_2) ; 
d4s  3 : factor (d4s_3) ; 


/*  DERIVE  V  FUNCTION  */ 

d_03  : matrix  ( [-d_0] ,  [mu_0] )  $ 

d_01  : matrix  (  [  d _ 0  ]  ,  [mu_0]  )  $ 

d~10:matrix( [-1] , [0] ) $ 


d_12 :matrix  ( [1] , [0] ) $ 
d_21 : matrix ( [-d_2] , [mu_2] ) $ 
d_23 : matrix { [d_2] , [mu_2] ) $ 
d_32:matrix([-l],  [0])$ 
d_30:matrix(  [1] ,  [0] )  $ 

r_10  rmatrix  ( [cos  (th_10)  ,  -sin  (th_10)  ] , 

[sin(th_10),  cos (th_10) ] ) $ 

3 _ 01 : - (3/8) * (transpose  <d_01) . (r_10 . d_10) )  + 

(1/8) * (transpose (d_03) . (r_10 . d_12) ) $ 
j_Cll :  rat  simp  ( j_01,  sin(th_10),  cos  (th_10)  )  $ 

r_20  imatrix  ( [cos  (th_20)  ,  -sin  (th_20)  ] , 

[sin(th_20),  cos (th_20) ] ) $ 
j_02:-(l/8) * (transpose (d_01)  . (r_20.d_23) ) - 

(1/8) * (transpose (d_03) . (r_20.d_21) ) $ 
j_02 : rat simp ( j_02,  sin(th_20),  cos(th_20))$ 

r_03  rmatrix  ( [cos  (th_03)  ,  -sin  (th_03)  ] , 

[sin(th_03),  cos (th_03) ] ) $ 

j _ 03 : — (3/8) * (transpose (d_30) . (r_03.d_03) ) + 

(1/8) * (transpose (d_32) . (r_03 . d_01 ) ) $ 
j_03 : ratsimp ( j_03,  sin(th_03),  cos(th_03))$ 

r_21  rmatrix  (  [cos  (th_21)  ,  -sin  (th_2l)  ]  , 

[sin(th_21),  cos (th_21) ] ) $ 
j_12 : - (3/8) * (transpose (d_12) . (r_21 . d_21) ) + 

(1/8) * (transpose (d_10) . (r_21 . d_23) ) $ 
j_12 : ratsimp (j_12,  sin(th_21),  cos(th_21))$ 

r_31 rmatrix  ( [cos  (th_31)  ,  -sin  (th_3l)  ] , 

[sin (th_31) ,  cos (th_31) ] ) $ 
j_13 : - (1/8) * (transpose (d_12) . (r_31 . d_30) ) - 

(1/8) * (transpose (d_10) . (r_3l . d_32) ) $ 
j_13 : ratsimp ( j _ 13 ,  sin(th_31),  cos(th_31))$ 

r_32  rmatrix ( [cos (th_32) ,  -sin (th_32) ] , 

[sin (th_32) ,  cos (th_32) ) ) $ 
j_23 : - (3/8) * (transpose (d_23) . (r_32 . d_32) ) + 

(1/8) * (transpose (d_21) . (r_32 .d_30) ) $ 
j_23 : ratsimp ( j_23,  sin(th_32),  cos(th_32))$ 

vrratsimp  ( j_01  + j_02+ j_03+ j_12+ j_13+ j_23,  mu_0)  $ 
v:ev(v,  th_03  -  -th_30,  th_32  -  th_30  -  th_20, 

th_31  -  th_30  -  th_10,  th_21  -  th_20  -  th_10) ; 

/**♦** **************************************/ 

/*  THE  TIRST  DERIVATIVE  OF  V  TUCTION  */ 

/*  W.R.T.  th_10  */ 

/*******************************************/ 
dlv: dif f (v,  th_10,  1)$ 

/*******************************************y 
/*  THE  FIRST  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  AT  */ 

/*  SYl^HETRIC  CONFIGURATION  */ 

/A******************************************/ 

dlvs:ev(dlv,  diff(th_20,  th_10,  l)—dls_2,  diff (th_30,  th_10, 
th_20-%pi,  fch_30— ths_10,  th_10»ths_l 0 ) $ 
dlvs :trigexpand (dlvs) ; 

/*******************************************/ 

/*  THE  SECOND  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  */ 

/**************************************** / 

d2v:diff (dlv,  th_10,  1)$ 


1) “dls  3, 


/*******************************************/ 

/*  THE  SECOND  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  AT  */ 

/*  SYMMETRIC  CONFIGURATION  */ 

ft******************************************/ 

d2va :ev (d2v,  diff(th_20,  th_10,  l)-dls_2,  diff(th_30,  th_10,  l)-dla_3, 
diff (th_20,  th_10,  2)-d2s_2,  diff(th_30,  th_10,  2)-d2s_3, 
th_20-%pi,  th_30— ths_10,  th_10-tha_10) $ 
d2vs rtrigexpand (d2va) ; 

/*******  ************************************/ 

/*  THE  THIRD  DERIVATIVE  OF  V  FUCTION  */ 

/*  M.R.T.  th_10  */ 

^** ***************************&*************  j 

d3v:diff (d2v,  th_10,  1) $ 

/♦A*****************************************/ 

/*  THE  THIRD  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  AT  */ 

/*  SYMMETRIC  CONFIGURATION  */ 

/*******************************************/ 

d3vs :ev (d3v,  diff(th_20,  th_I0,  l)-dls_2,  diff<th_30,  th_10,  l)-dls_3, 
diff (th_20,  th_10,  2)«d2s_2,  diff(th_30,  th_10,  2)-d2a_3, 
diff (th_20,  th_10,  3)-d3s_2,  diff(th_30,  th_10,  3)-d3a_3, 
th_20«%pi,  th_30— tha_10,  th_10-ths__10)  $ 
d3vs:trigexpand (d3vs) ; 

/*******************************************/ 

/*  THE  FOURTH  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  */ 

/A******************************************/ 

d4v : diff (d3v,  th_10,  1)$ 

/*******************************************/ 

/*  THE  FOURTH  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  AT  */ 

/*  SYMMETRIC  CONFIGURATION  */ 

/*******************************************/ 

d4va:ev(d4v,  diff (th_20,  th_10,  1)— dla_2,  diff(th_30,  th_10,  l)—dla_3, 
diff (th_20,  t h_l 0 ,  2)-d2a_2,  diff(th_30,  th_10,  2)-d2a_3, 
diff  (th_20,  th_10,  3)«d3a_2,  diff(th_30,  th_10,  3)-d3a_3, 
diff (th_20,  th_10,  4)-d4s_2,  diff(th_30,  th_10,  4)-d4s_3, 
th_20-%pi,  th_30—tha_10,  th_10— tha_10) $ 
d4va :trigexpand (d4vs) ; 


/*******************************************/ 

/* 

THE  SECOND  DERIVATIVE  OF 

V  FUCTION 

V 

/* 

W.R.T.  th_10  AND 

*/ 

/* 

THE  FIRST  DERIVATIVE  OF  V 

FUNCTION 

*/ 

/* 

W.R.T.  mu  0 

*/ 

/*******************************************/ 
d3vmu:diff  <d2v,  mu_0,  1)$ 

/*******************************************/ 

/*  THE  SECOND  DERIVATIVE  OF  V  FUCTION  */ 

/*  W.R.T.  th_10  AND  */ 

/*  THE  FIRST  DERIVATIVE  OF  V  FUNCTION  */ 

/*  W.R.T.  mu_0  */ 

/*  EVALUATED  AT  SYMMETRIC  CONFIGURATION  */ 

/*******************************************  j 

d3vtnua : ev  (d3 vmu,  diff  (th_20,  th_10,  l)-dla_2,  diff(th_30,  th_10,  l)-dla_3 
diff (th_20,  th_10,  2)-d2a_2,  diff(th_30,  th_10,  2)-d2a_3 
th_20-%pi,  th_30— tha_10,  th_10-tha_10 )  $ 
d3vmua : trigexpand (d3vmua ) ; 


/a***********************************************/ 

/*  FIND  mus_0  AT  WHERE  THE  BIFURCATION  HAPPENS  */ 
/************************************************/ 
linsolve ( [d2vs] ,  [mu_0] ) ,globalsolve:true$ 
mus_0:trigexpand (mu_0) ; 

/*******************************************************/ 


/*  AT  BIFURCATION  POINTS  THE  FOLLOWING  EQUATIONS  */ 
/*  SHOULD  HOLD  */ 
/*  g(ths_10,  mus_0)  “0  */ 
/*  Dg_th_10  (ths_10,  mus_0)  “0  */ 
/*  DDg_(th_10,  th_10) (ths_10,  mus_0)  —  0  */ 
/*  WHERE  g-Dv/Dth_10,  ths_10  IS  th_10  AT  */ 
/*  SYMMETRIC  CONFIGURATION  */ 


/a******************************************************/ 
trigsimp (ev (dlvs,  mu_0— mus_0) ) ; 
trigsimp (ev (d2vs,  mu_0— mus_0) ) ; 
trigsimp (ev (d3vs,  mu_0-mus_0) ) ; 

/*+ ★★* **★*******★*★***★*****★*★*★★**★★**★***★★*★★★**★*★★/ 

/*  AT  BIFURCATION  POINTS  THE  FOLLOWING  EQUATIONS  */ 

/*  SHOULD  HOLD  */ 

/*  DDDg_(th_10,  th_10,  th_10) (ths_10,  mus_0)  not-  0  */ 

/*  DDg_(th_10,  mu_0) (ths_10,  mus_0)  not-  0  */ 

/♦a*****************************************************/ 

d4vs  :ti:igexpand  (trigsimp  (ev  (d4vs,  mu__0— mus_0)  )  )  ; 
d3vmus : trigexpand (trigsimp (ev (d3vmus,  mu_0— mus_0)  )  )  ; 

/*******************************************************/ 

/*  AT  LEADING  FORM  */ 

/*  MAKE  d4vs  AND  d3vmus  AS  TUCTIONS  OF  d_0  AND  d_2  */ 

/*******************************************************/ 
musl_0 : ev (mus_0,  sin (ths_10) — sqrt (4 . - (d_2-d_0) **2) /2, 
cos  (ths_10)  —  (d_0-d_2)  /2)  $ 
musl_0 : ratsimp (musl_0) ; 

d4vslead :  ratsimp  (ev  (d4vs,  sin  (ths_10)  —sqrt  (4  .  -  (d_2-d_0)  **2)  /2, 

cos (ths_10) — <d_0-d_2) /2, 
mus_0— musl_0) ) ; 

d3vmuslead:  ratsimp  (ev  (d3vmus,  sin  (ths_10)  —sqrt  (4  .  -  (d_2-d_0)  **2)  12, 

cos  (ths_10)  —  (d_0-d_2)  /2, 
mus_0-musl_0) ) ; 

/*******************************************************/ 

/*  AT  LAGGING  FORM  */ 

/*  MAKE  d4vs  AND  d3vmus  AS  FUCTIONS  OF  d_0  AND  d_2  */ 

/*********★*********************************************/ 
mus2_0 : ev (mus_0,  sin (ths_10) -+sqrt (4 . - (d_2-d_0) **2) 12, 
cos (ths_10) — (d_0-d_2) / 2) $ 
mus2_0 : ratsimp (mus2_0) ; 

d4vslag:  ratsimp  (ev  (d4vs,  sin  (ths_10)  — +sqrt  (4  .  -  (d_2-d_0)  **2)  /2, 

cos  (ths_10)  —  (d_0-d_2)  /2, 
mus_0-mus2_0) ) ; 

d3vmusl.ag:  ratsimp  (ev  (d3vmus,  sin  (ths_10 ) — +sqrt  (4  (d_2-d_0)  **2 )  12, 

cos  (ths_10)  —  (d_0-d_2)  12, 
mus  0^mus2  0)  )  ; 


/*  AN  EXAMPLE  */ 

/********************* / 

ev(musl_0,  d_0-4/3,  d_2— 1); 
ev(mus2_0,  d_0— 4/3,  d_2— 1); 

ev(d4vslead,  d_0— 4/3,  d  2—1); 
ev (d3vmuslead,  d_0-4/3,  d  2-1); 


ev(d4vslag,  d_0-4/3,  d_2-l) ; 
ev (d3vmuslag,  d_0— 4/3,  d_2—l) 

apply (closefile, [bifur_out] ) ; 


